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Abstract 

Accurate knowledge of cross sections for pion production in proton-proton col- 
lisions finds wide application in particle physics, astrophysics, cosmic ray physics 
and space radiation problems, especially in situations where an incident proton is 
transported through some medium, and one requires knowledge of the output par- 
ticle spectrum given the input spectrum. In such cases accurate parametrizations 
of the cross sections are desired. In this paper we review much of the experimen- 
tal data and compare to a wide variety of different cross section parametrizations. 
In so doing, we provide parametrizations of neutral and charged pion cross sec- 
tions which provide a very accurate description of the experimental data. Lorentz 
invariant differential cross sections, spectral distributions and total cross section 
parametrizations are presented. 
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1. Introduction 

Pion production in proton-proton collisions has been extensively studied over many 
years, and has now reached the point where this knowledge finds useful applications in a 
variety of detailed below. 

1. Two important types of particle detectors are the hadronic and electromagnetic 
calorimeters Jl|, in which an electromagnetic or hadronic shower is initiated by a high 
energy incoming particle. From a Monte-Carlo simulation of the shower, one is able to 
deduce important characteristics of the incoming particle such as its energy and identity. 

2. The primary cosmic rays can be detected by a variety of methods, depending on 
the incident energy. For the very high energy cosmic rays, where the flux is relatively low, 
the extensive air showers (EAS) 0, |3|, f|] provide the most convenient means of detection. 
The EAS is analogous to the hadronic or electromagnetic calorimeter used in particle 
physics, but with the Earth's atmosphere being the active volume in which the shower 
develops. The EAS has both electromagnetic and hadronic components, and similar to 
the calorimeter, the energy and identity of primary cosmic ray nuclei can be deduced via 
Monte-Carlo simulation of the showers @, |3|. 

3. In long duration human space flights, such as a mission to Mars, the radiation 
levels induced by galactic cosmic rays can exceed exposure limits set for astronauts 0, ||. 
In determining the radiation environment inside a spacecraft one needs to transport the 
exterior cosmic ray spectrum through the spacecraft wall in order to determine the interior 
radiation spectrum. 

4. In gamma ray |7|, |8|] and high energy neutrino astronomy |9|, [TI| , the diffuse back- 
ground radiation is due in large part to the gamma rays and neutrinos produced in proton 
collisions with the protons in the interstellar medium. In addition, pion production from 
proton-proton collisions finds applications in the calculation of gamma ray emission from 
the accretion disk around a black hole jnj . 

In all of the above applications it is crucial to have an accurate knowledge of the 
cross sections for pion production in proton-proton collisions. In addition, most of the 
applications mentioned above require solving the transport equations which determine 
the particle spectrum on one side of a material (active volume of calorimeter, Earth's 
atmosphere, spacecraft wall or interstellar medium) given the incident particle spectrum. 
Use of pion production cross sections in such transport codes requires that the cross 
section be written in a simple form. The transport codes have many iterative loops, 
which will take too much computer time if the cross section formulae also contain many 
iterative loops. Thus it is most advantageous if one can write down simple formulae which 
parameterize all of the experimental data on pion production cross sections. That is the 
aim of the present work. 

In this paper we provide simple algebraic parametrizations of charged and neutral pion 
production cross sections valid over a range of energies. The cross sections we provide are 
Lorentz- invariant differential cross sections (LIDCS), lab frame spectral distributions (i.e. 
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energy differential cross sections) and total cross sections, because these are the types 
of cross sections most widely used in transport equations. Many authors have presented 
such parametrizations before, but the problem is to decide which authors are correct and 
whether a particular parametrization applies only to a limited data set or is valid over a 
wider range. In the present work, we have performed an exhaustive data search and have 
compared as many different parametrizations as possible to as much data as possible, so 
as to provide definitive conclusions as to which is the most accurate parametrization to 
use. All of this is discussed more extensively below. 

The cross sections discussed in this paper are for inclusive pion production in proton- 
proton collisions, i.e. the reactions considered are p + p — > n + X, where p represents a 
proton, 7i represents a pion, and X represents any combination of particles. An extensive 
search for LIDCS data was performed, and the data was used to compare all available 
parametrizations. An extensive set of data was used in these comparisons, but only a 
few data points are graphed in this paper due to space considerations. A method for 
generating parametrizations for these cross sections is also described and applied to n° 
production. Spectral distribution and total cross section formulae were not developed di- 
rectly because of a lack of data. Instead, the most successful LIDCS parametrizations were 
first transformed into lab frame spectral distributions by numerical integration. These 
spectral distributions were parameterized and then numerically integrated to generate 
lab frame total cross sections. Finally, the total cross sections were compared to available 
data and parameterized as well. This procedure is discussed, and the parametrizations of 
the numerical results are given. Multiple checks of the accuracy of all results were made, 
and some of these are presented. 

Notation: Starred quantities (e.g. 9*) refer to the quantities in the center of mass 
(COM) frame, while unstarred quantities (e.g. 9) refer to the quantities in the lab frame. 
= Lorentz-Invariant Differential Cross Section (LIDCS) 

|§ = spectral distribution = 2irp J° max dBE^wnO 

a = total cross section = 2vr f n max d9 f Pmax dpE^f f sin9 = \l max %dE 

P p is the proton momentum. 
m p is the proton mass. 
m n is the pion mass. 

\fs is the magnitude of the total four momentum, and is equal to the total energy in the 
COM frame. 

Ti a b is the lab frame kinetic energy of the incoming proton. 
T is the pion kinetic energy. 
E is the pion total energy. 

9 is the angle of pion scattering with respect to the direction of the incident particle. 
p is the pion momentum. 

P± = Pt is the pion transverse momentum (p± = psin9). 

Pmax is the maximum possible momentum the scattered pion can have for a given y/s. 
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2. Comparison of Lorentz Invariant Differential Cross Sections 



The object has been to determine an accurate parameterization for inclusive LIDCS, 
which can be confidently applied to regions where no experimental data is available. For 
example, the parametric equation would need to be extrapolated to energies lower than 
those for which data are available, if the formulae were to be used for the purpose of 
developing radiation shielding materials. The most convenient formulae are those that 
are in closed form, since they are easily used, and take relatively little CPU time in 
numerical calculations. Some of the formulae that were considered as representations 
of the LIDCS were not in closed form, but included tabulated functions of energy (i.e. 
numerical values were given for specific energy values rather than a functional form). 
When comparing parameterizations, closed form expressions were given precedence over 
other equally accurate formulae. 

The invariant single-particle distribution is defined by 

^ ^ c d 3 p c d 3 p p 2 dpdVL ^ ^ 

where Jj^- is the differential cross-section (i.e. the probability per unit incident flux) 
for detecting a particle C within the phase-space volume element d 3 p c . A and B are the 
initial colliding particles, C is the produced particle of interest, and X represents all other 
particles produced in the collision. E is the total energy of the produced particle C, and 
Q is the solid angle. This form is favored since the quantity is invariant under Lorentz 
transformations. 

The data for pion production in proton-proton interactions is primarily reported in 
terms of the kinematic variables 8*, ^/s, p±, which are respectively, the center of mass 
(COM) frame scattering angle of the pion, the invariant mass of the entire system, and 
the transverse momentum of the produced pion. a/s is a Lorentz invariant quantity, and 
is equal to the total energy in the COM frame. p± = p*sin#*, where p* is the COM 
momentum. p± is invariant under the transformation from the lab frame to the COM 



frame. (See [|2| for a more detailed discussion of kinematic variables.) In the following 



discussions, all momenta, energies, and masses are in units of GeV. 
2.1 Neutral Pions 



Busser et al. |T3j have fitted the LIDCS data obtained in the reaction p + p — > 7r° + X, 
where p represents a proton, 7r° represents the neutral pion produced, and X represents 
all other produced particles, to an equation of the form 

with A = 1.54 x 10~ 26 , n = 8.24, and b = 26.1. This equation is based on a specific set of 
experimental data with all measurements taken at 6* ~ 90°, and was originally intended 
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only for pions with high p±. Comparison of this parameterization with data available 
from other experiments [14j] - indicates that the global behavior of the invariant cross 



section cannot be represented by a function of this form. See Figures [L] - |3] for some 



examples of data. The parameterization of Busser et al. |L3| was not plotted because the 
cross section is much too small compared to the data in the p± ranges covered by the 
graphs. 



The following form has been used by Albrecht et al. [^] to represent neutral pion 
production. 

^ = c{^-r (3) 

d 3 p V+PO 

where C, n, and po are free parameters. This equation only has dependence on p± where 
as the data [H - ||20|| , some of which is shown in Figures |I] - also has dependence on 



y/~s and 6*. This form is therefore not general enough to represent all the data. 



Ellis p2[ have favored a representation for the invariant cross section of the form 



E^ = A(pl + M*r N/2 f(x±,0*) (4) 

where f(x±, 6*) = (1 — x±) F , N and F are free parameters and the scaling variable x± is 
given by x± = ~ p* max = [{s + m 2 — 4m 2 ) 2 /4s — m 2 ] 2 , where and m p are the 
mass of the neutral pion and the proton respectively, is the maximum pion momentum 
allowed. The outline of this basic form has been used by Carey et al. in fitting the 
invariant cross section for the inclusive reaction p + p —* tc° + X Their representation 
is given by 

E^ = A(pl + 0.86r^(l-x* R y (5) 

where x* R = -f — is the radial scaling variable and the normalization constant A has 
been determined as A ~ 5. This parameterization accurately reproduces the data for 
measurements taken at 6* = 90° and ^fs > 9.8 GeV , but does not agree well with the 
data for lower energies (yfs = 7 GeV). The disagreement at lower energies can be seen in 
Figure [I]. 

Another problem with this parameterization becomes apparent, when one considers 
that integration over all allowed angles and outgoing particle momenta should yield the 
total inclusive cross section. The details of this calculation appear in Section 3. A 
comparison of the experimentally determined total cross section data from Whitmore 
with the results of the numerical integration of equation shows that the total cross 
section is greatly underestimated by Carey. See Figure f§. 



Stephens and Badhwar |20| obtained data from the photon cross sections given by 
Fidecaro [fLUR . The Fidecaro data was taken at incident proton kinetic energy T/ a b = 23 
GeV and p± = 0.1 GeV - 1.0 GeV. (Note: No error was listed by Fidecaro et al. [14] for 



pion production. Error bars of ten percent were added to the data on the figures, since this 
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level of error was standard for most of the other data. Also, Stephens uses the notation 
E p instead of TJ afe .) Figures p] - [3], show common examples of the accuracy of Badhwar 
and Stephens' parameterization's fit to the data. The following is the parameterization 
of the 7r° invariant cross section proposed by Stephens and Badhwar |^0| ; 

E^- = Af(T lab )(l-xyexp(-BpJ(l+4m 2 p /s)) (6) 

where 

x 

q 

f(Tlab) 

and 

A = 140,5 = 5.43, Ci = 6.1, C 2 = 3.3, C 3 = 0.6 with xt = and pi = p* cos6*. 

The Stephens-Badhwar parameterization was found to be the best of the previously 
listed representations, because it accurately reproduces the data in the low p± region, 
where the cross section is greatest (see Figures [I] - |3|), and its integration yields accurate 
values for the total cross section (see Figure 0). This equation is, however, a poor tool for 
predicting values of the invariant cross section for p± > 3 GeV, where the value predicted 
underestimates experimental data by up to ~ 10 orders of magnitude (see Figures || and 

D- 

No parameterization currently exists that accurately fits the global behavior of the 
LIDCS data. Previous equations have suffered from being too specific to a particular set 
of experimental data, or from failing to reproduce the total cross section upon integration. 
It is for these reasons that a new parameterization is desired, one that correctly predicts 
all available data while maintaining the essential quality of correctly producing the total 
cross section upon integration. 

The approach that has been adopted in the present work is to assume the following 
form for the invariant cross section 

E^- = (sin e*) D{ ^ ,P± ' et) F(^s, p±,9* = 90°) (7) 
d A p 

The motivation for an equation of this form is that as the angle decreases, the cross section 
decreases very slowly at lower p± values. The approximation that was made in deriving 
the above equation is that as p± — > 0, the cross section is assumed to be independent of 
angle. 



C 1 -C 2P± + C 3 p* ± 
yjl + Am 2 p /s 
= (l + 23T^- 6 )(l-4m;/ S ) 2 
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Under the assumption that the invariant cross section can be fitted by equation (|7]), 
the program goes as follows. Find a representation for the cross section as a function 
of energy yfs and transverse momentum p± from experimental data taken at 9* = 90°. 
F(- s /s,p±) is then completely determined, because (sin#*) D is unity at 9* = 90°. 

At 9* = 90° the data is well represented by 

E^-(9* = 90°) = F(V~s,P±), 

with 

F{V5,p ± ) = hi(-j£-)G(q,p ± ), (8) 

where q = and the COM pion production threshold energy y/s min = 2m p + m^. The 
function, 

E€£(9* = 90°) 

(•■[<!■ l>.) 



ln( r < 

V min 

was parameterized as 

G{q, p±) = exp{fci + k 2 p± + hq' 1 + k A p\ + k 5 q~ 2 + k & p±q~ l + k 7 p\ + k 8 q~ 3 

+k 9 p±q' 2 + hoplq' 1 + k u pl 3 } (9) 

with kt = 3.24, k 2 = -6.046, k 3 = 4.35, k A = 0.883, k 5 = -4.08, k 6 = -3.05, k 7 = 
-0.0347, k 8 = 3.046, k 9 = 4.098, k 1Q = -1.152, and k n = -0.0005. The parameters 
k% — k w were obtained using the numerical curve fitting software Table Curve 3D v3 [P^j 
and the eleventh term was added to modify the low p± behavior of the parameterization. 

With F(y/s,p±) determined, the function D(y/s,p±,6*) is the only remaining un- 
known. Solving for D yields 

D ^ 8 ^ e ) = ( 10 ) 

Equations (§) and (||) were then used in equation (" pT0|) to calculate values of D(y/s,p±, 9* 
If the function D is independent of angle, then equation flip]) could be determined for any 
fixed angle, 9* ^ 90°. Data were compared for a range of angular values, and this data 
revealed that the function D is not independent of angle. The angular dependence turned 
out to be of the form (sin#*) _a45 , and 

D(V~s,P±,9*) = (sm9T 0A5 [ciP c l(V~s) C3 + c^ + ^ + —] (11) 

ys y s s 
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with ci = 205.7, c 2 = 3.308, c 3 = -2.875, c 4 = 10.43, and c 5 = 0.8. 

The final form of our resultant parameterization for the neutral pion invariant cross sec- 
tion in proton-proton collisions is equation (Qj with F(p±, y/s) given in equation (^), 
G(q,p±) given in equation^), and D{p A 



's,9*) given in equation (11). This form is 
accurate over a much greater range of transverse momentum values than those covered by 
previous representations. Figures [I] - |3] show a few comparisons. A much more extensive 
set of data was used in the development and comparison of the parameterizations, but 
they are not shown in this paper due to space considerations. For the low transverse 
momentum region where the cross section is the greatest, the fit is quite similar to that 
of Stephens et al. |2(J . Also, Figure £| shows that both formulae (|6| and 0) integrate to 



approximately the same total cross section, which is in agreement with the data from 
Whitmore et al. [FJ|]. (Equation [7] integrated into a total cross section is denoted as 
Kruger in Figure A more complete comparison of the integrated total cross section to 
data is given by Stephens et al. pD| . Note however that equation (|7]) was based mainly 
on the data from |TJ| - |20[] . Equation (|^) could therefore give unpredictable results in 
regions not included in those data sets, particularly for very low transverse momentum 
or yfs > 63 GeV. 



8 



2.2 Charged Pions 

The available data for charged pions, is less extensive than ir° data. There is therefore 
a higher degree of uncertainty in LIDCS for charged pions. Integration of a LIDCS to get a 
total cross section and comparison of the results to total cross section data, allows a check 
of the global fit of a parametrization. This check was made for charged as well as neutral 
pions, but due to a lack of data, it is more important for charged pions. Parametrizations 
that do not integrate to the correct total cross section can be ruled out, even if the LIDCS 
data is well represented, because the global behavior of the parametrization cannot be 
accurate. However, producing a correct total cross section upon integration does not 
necessarily imply that the global behavior of the parameterization is correct. A tighter 
constraint could be placed on possible LIDCS parametrizations, if more measurements 
were made. If the spectral distribution is measured at three different values of pion energy 
for two different proton collision energies, the general behavior of the spectral distribution 
could be checked. The angular dependence of LIDCS parametrization could then be tested 
by integrating over angle, and comparing the results to the spectral distribution data. For 
the purposes of space radiation shielding, measurements at proton lab kinetic energies of 
3 GeV and 6 GeV, and pion lab kinetic energies of 0.01 GeV, 0.1 GeV, and 1 GeV would 
be useful, because this is the region with both a large cross section, and large Galactic 
Cosmic Ray fluxes. With these facts in mind, a comparison of LIDCS parametrizations 
with data from Jl6|, |25], [27], for charged pion production follows. 



A parametrization for n of the form 



E j^ = Ae M-B P l) (12) 



has been given by Albrow et al. f29j where A and B are tabulated functions of x R = -f — . 
A and B are given only for x R = 0.18, x R = 0.21, and x R = 0.25 which limits the 
usefulness of this parametrization. 

Alper et al. [55" have fitted the data for both ir + and 7r~ production to the following 
form 

E^- = Aexp(-B P± + Cpl) exp(-ZV) (13) 

where y is the longitudinal rapidity, and A, B, C, and D are tabulated functions of s, that 
are also dependent on the type of produced particle (n + or (Note that at 6* = 90° 
we have y — 0.) The fit to the data is excellent for low transverse momentum, as can be 
seen in Figures |8| and [TT], but these figures also show that this form has an increasing cross 
section for high p±, which contradicts the trend in the data. Also, there are different sets 
of constants for each different energy, which makes a generalization to arbitrary energies 
difficult. 



Parametrizations done by Carey et al. |30| and Ellis et al. [22] have a similar form 



although Carey's was applied only to n . Both underestimate LIDCS for low p±, where 
the cross section is the largest (see Figures 0-0). 
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The following is Carey's parametrization 



d 3 rr 

E— (tt-) = N(pl + 0.86)" 4 - 5 (l - ^) 4 (14) 
where N=13 is the overall normalization constant, and x* R = -§ — ~ 2 j=. 

Pmax V s 

The following is Ellis's parametrization which was applied to both 7r + and tt~ production 
at 6* = 90°. 

E^- = A( P l + M 2 )- N/2 (l-x±f (15) 

where M, N, F are given constants. A is an unspecified overall normalization for which 
we used A = 13, and x± = ^p— ps 

The most successful LID(?S parametrization available for charged pion production was 



found to be the one developed by Badhwar et al. [31 



d 3 a = A(l-x) q r-Bp, /(i +4m 2/ s )i n6] 
d 3 p (l+4m 2 p /sY 1 } 

where q is a function of p± and s, such that 

9 = (C 1 + C 2 p ± + C 3 p 2 ± )/(l + 4m 2 p /s) 1 / 2 

and 

x « [xf + -{p 2 ± + m 2 n )}i 
" s 

Here x], = m 2 A. For vr + , A = 153, 5 = 5.55, C l = 5.3667, C 2 = -3.5, C 3 = 0.8334, 

II Pmax V s 

and r = 1. For vr", A = 127, 5 = 5.3, <7 a = 7.0334, C 2 = -4.5, C 3 = 1.667, and r = 3. 
This form is accurate for low transverse momentum (Figures |7HT2"D, which is the most 
important region for radiation shielding due to the large cross section. It is also in closed 
form, so that extra numerical complexities do not have to be considered. A comparison to 
a few data points, shown in Figure [H^, demonstrates that it integrates to the correct total 
cross section. A more detailed comparison of the integrated cross section to experimental 



data is given by Badhwar et al. Because of its relative accuracy and simplicity, this 
parametrization was integrated to get total cross sections and spectral distributions for 
charged pions. 



Mokhov et al. [32] have also developed the following formulae for both tt + and 7r 
production. 

E^f = A(l - -^) B eM-^rWMVM (17) 
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where 



(1 -D)exp( 

0.2625 
{p 2 ± + 0.87) 4 



+ Dexp(-Fp5. 
for p ± > 0.933 GeV 



for Pl < 0.933 GeV 



and 



V 2 = 0.7363 exp(0.875p ± ) for p ± < 0.35 GeV 
= 1 for p ± > 0.35 GeV 

with A=60.1, B=1.9, and C=0.18 for tt+, A=51.2, B=2.6, and C=0.17 for vr", and D=0.3, 
E=12, and F=2.7 for both ir + and ix~ . Figures show that the formula of Badhwar 
has a better fit to the data in the low p± region where the cross section is the largest. 
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3. Spectral Distributions and Total Cross Sections 



3.1 Method of Generating Other Cross Sections from a LIDCS 

While LIDCS contain all the necessary information for a particular process, sometimes 
other cross sections are needed. For example, one dimensional radiation transport requires 
probability density distributions that are integrated over solid angle. These quantities are 
calculated in terms of spectral distributions and total cross sections rather than LIDCS, 
but with accurate parametrizations of LIDCS, formulae for both spectral distributions 
and total cross sections can be developed. LIDCS for inclusive pion production in proton- 
proton collisions contain dependence on the energy of the colliding protons (a/5), on the 
energy of the produced pion (T n ), and on the scattering angle of the pion (9). Total cross 
sections a, which depend only on y/s, and spectral distributions j^, which depend on yfs 
and T n can be extracted from a LIDCS by integration. If azimuthal symmetry is assumed, 
these cross sections take the following forms 

% = 2np dOE^smd (18) 

dE J d A p 

a = 2tt / d6 dpE— - r 19 



m i n 



where 9 max , p m ax, and p min are the extrema of the scattering angle and momentum of the 
pion respectively, and m n is the rest mass of the pion. 

In the Center of Mass (COM) frame these extrema can easily be determined. Using 
conservation of momentum and energy, one can easily show that 

p2= (s + m l-s^ _ ml (20) 

where s x is the square of the invariant mass of the sum of all particles excluding the pion, 
and p is the magnitude of the three-momentum of the pion. The independence of p on 
9 implies that 9 can take on all possible values (ie. 9 max = ir), and the symmetry of the 
COM frame implies that p m in = 0. For a given value of s, it is obvious that momentum 
is a maximum when s x is a minimum. An invariant mass is a minimum, when it is equal 
to the square of the sum of the rest masses of the particles in question. Momentum is, 
therefore, a maximum when s x is the square of the sum of the least massive combination 
of particles that can be produced while still satisfying all relevant conservation laws. For 
the reaction p + p — > ir + x, we have s x ~ 4mp, where a subscript p represents a proton. 

If a Lorentz transformation is applied to the maximum COM momentum, the integra- 
tion limits can be determined in other frames. Byckling and Kajantie have shown that 
by transforming to the lab frame, the following formula can be obtained \Y1\ 



p± = \p a E* max yG C os9±(E a + m p )Jsp% ax -mlplsm 2 9}[s+plsm\9)]- 1 (21) 
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where starred quantities are COM variables, and unstarred quantities are either lab or 
invariant variables, m p is the rest mass of a proton; p a is the magnitude of the momentum 
of the projectile proton, and p + = p max is the maximum pion momentum. The greater of 
the two quantities p~ = p m in and is the minimum pion momentum, and the maximum 
scattering angle can be determined by the requirement that p ± be real. This requirement 
implies that the quantity under the square root must be greater than or equal to 0. Solving 
for 9 max then gives the formula 



With the limits of integration determined, a LIDCS can be turned into a total cross section 
or a spectral distribution by numerical integration. This procedure will, however, give 
discrete "data" points; not closed form expressions. Parametrizations of this numerical 
data are needed, if relatively simple formulae for these cross sections are desired. This 
process was completed for all three pion species, and the corresponding formulae are 
listed in the next section. It should be noted that the accuracy of these parametrizations 
is limited to that of the original LIDCS. 




p a m. 



(22) 



'7T 
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3.2 Parametrizations 



The surface parametrizations for the spectral distribution as a function of incident pro- 
ton kinetic energy in the lab frame {Ti a b) and the lab kinetic energy of the produced pion 
(T n ) has been completed by numerically integrating LIDCS charged pion parametrizations 
due to Badhwar et al. (equation |HJ and the neutral pion cross section both from 
Stephens and Badhwar (equation ^) P0| , and from equation [7]. The numerical integration 
routines were checked by computing total cross sections in both the lab and COM frames 
and comparing the results. Since total cross section is invariant under the transformation 
between these two frames, the results should be the same in both frames. In order to 
accurately fit the integration points for low energies, it has been necessary to consider 
two regions of the surface and to determine representations for them individually. For 
each of the three pions, the two regions consist of laboratory kinetic energies (Tj a &) from 
0.3 GeV to 2 GeV and from 2 GeV to 50 GeV. Using the following parametrizations in 
energy regions other than the region listed above could give unpredictable results since 
the formulae were not tested there. 

The neutral pion spectral distribution for the range 0.3 - 2 GeV is represented by the 
following equations: 



F 2 



A X T£* + A 3 T A4 



lab 



exp(v4 5 + 



A f 



+ A 7 T% + A 9 T A ™ + AuT. 



lab 



[%)iab = {A u ^+A u ew(A w VTZ + A 17 T^Tfr))Tf« 



(23) 



with constants Ai given in Table 1. 



Ai = 6.78 x 1CT 1U 


As = -1.75 


Am = 0.25 


A 2 = -2.86 


A 9 = -32.1 


Aw = -39.4 


A 3 = 1.82 x 1CT S 


Aio = 0.0938 


An = 2.88 


A4, = -1.92 


An = -23.7 


Ais = 0.025 


A 5 = 22.3 


A12 = 0.0313 


Aw = 0.75 


A 6 = 0.226 


Ai 3 = 2.5 x 10° 




A T = -0.33 


A 14 = 1.38 





Table 1: Constants for equation (23) 



Bi = 1.3 x lO" 1 " 


B 8 = -1.25 


Bis = 60322 


B 2 = -2.86 


Bg = -33.2 


Bl6 = 1.07 


B 3 = 4.27 x 10~ y 


B 10 = 0.0938 


B 17 = -67.5 


B 4 = -2.4 


Bn = -23.6 




B 5 = 22.3 


B 12 = 0.0313 




B 6 = -1.87 


B13 = 2.5 x 10" 




B 7 = 1.28 


B 14 = 0.25 





Table 2: Constants for equation (24) 



The neutral pion spectral distribution for the range 2-50 GeV is represented by the 
following equations: 



B X T** + B 3 T t B J 
B 6 



F\ = exp(B 5 + 



lab 



[^) lab = B 13 T^^ + B 15 T^exp(B 17 VT w ) 



(24) 
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with constants given in Table 2. 

The positively charged pion spectral distribution for the range 0.3 - 2 GeV is repre- 
sented by the following equations: 



, do 



F 2 
Fx 

)lab 



exp(C 5 + -Jg- + C 7 T% + C 9 T^ + CnT^T^ + C 14 In T lab ) 



lab 



= C 15 T?" ^ + C 17 T?" exp(C 19 + C 20 y/lQ 



with constants given in Table 3. 



(25) 



Ci = 2.2 x 10~ 8 


C 8 = -1.75 


Cib = 2.5 x 10" 


C 2 = -2.7 


C 9 = -29.4 


C*i6 = 0.25 


C 3 = 4.22 x 10" ' 


Cio = 0.0938 


C*i7 = 976 


C 4 = -1.88 


Cn = -24.4 


Cis = 2.3 


C 5 = 22.3 


C12 = 0.0312 


C*i9 = -46 


C 6 = 1.98 


C13 = 0.0389 


C20 = -0.989 


CV = -0.28 


C*i4 = 1.78 





Table 3: Constants for equation (25) 



D 1 = 4.5 x 10" 11 


D 7 = -35.3 


D13 = 60322 


D 2 = -2.98 


D 8 = 0.0938 


D14 = 1.18 


D 3 = 1.18 x 10~ 9 


A) = -22.5 


D15 = -72.2 


D 4 = -2.55 


D10 = 0.0313 


Die = 0.941 


D 8 = 22.3 


D n = 2.5 x 10" 


D 17 = 0.1 


D 6 = -0.765 


D 12 = 0.25 





Table 4: Constants for equation (26) 



The positively charged pion spectral distribution for the range 2-50 GeV is represented 
by the following equations: 

F 2 = D 1 T?*+D 3 T l % 

F 1 = exp(D 5 + -^ + D 7 T°«+D 9 T^») 

Vi lab 

{^)iab = D n T^^ + D n T^eMDi 5 VT w + D w T^) (26) 

with constants Di given in Table 4. 

The negatively charged pion spectral distribution for the range 0.3 - 2 GeV is repre- 
sented by the following equations: 

F 2 = dTp + GsTft 

F x = ex P (G 5 + ^^ + G 7 T^ + G 9 T^) 

vTi ab 

{%U = T^(G 12 ^ + G 13 eMGuVT n )) (27) 



with constants Gi given in Table 5. 
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Gi = 1.06 x 10- y 


G 6 = -1.5 


Gn = 0.25 


G 2 = -2.8 


G 7 = -30.5 


G12 = 2.5 x 10 b 


G 3 = 3.7 x 1CT H 


Gs = 0.0938 


G13 = 7.96 


G 4 = -1.89 


G 9 = -24.6 


G14 = -49.5 


G 5 = 22.3 


Gio = 0.0313 





Table 5: Constants for equation (27) 



Hi = 2.39 x 10- 1,J 


H 7 = -31.3 


H 13 = 60322 


H 2 = -2.8 


ifg = 0.0938 


H M = 1.1 


H 3 = 1.14 x 10- s 


H 9 = -24.9 


H 15 = -65.9 


Hi = -2.3 


Hio = 0.0313 


H\e = -9-39 


H 5 = 22.3 


fill = 2.5 x 10 b 


iTi7 = -1.25 


H 6 = -2.23 


H12 = 0.25 





Table 6: Constants for equation (28) 



The negatively charged pion spectral distribution for the range 2-50 GeV is repre- 
sented by the following equations: 

F 2 = H X T?* + H*T5Z 

F x = e X p{H s + -j^-+H 7 T* + H 9 T* a ) 

VI lab 

{%U = HnT^ + H^ex^^n + ^eT^) (28) 

with constants Hi given in Table 6. 

Total inclusive cross sections are represented by the following equations. 



a* = (0.007 + 0.1^2^ + ^)- 1 (29) 

1 lab J- i a b 



ln(T lab ) 0.162, , , , 

a w+ = (0.00717 + 0.0652 K labJ + -^Y 1 (30) 

lab J- lab 

0.0846 0.577 N , , , 

a ^ = (0.00456 + — — + —^y 1 (31) 

1 lab 1 lab 

For neutral pions, spectral distributions and total cross sections that were based on 
our own parametrization given in equation (Qj were also developed. The formula for the 
spectral distribution was not divided into two regions, and it is much simpler than the 
previous formulae. 

where K x = -5.8, K 2 = -1.82, K 3 = 13.5, K A = -4.5. 

Because equation (|7|) and Stephens LIDCS parametrization integrate to nearly the 
same total cross section (see Figure f|), separate total cross section parametrizations are 
not necessary (i.e. use equation |29|) . 

3.3 Discussion of Figures 

As discussed previously, Figures [T] - [3] show a comparison of LIDCS parametrizations 
for 7T° production of Carey et al. (equation ||) |1| , Stephens et al. (equation ||) |2(J , and 
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equation (|7]) plotted with data from jTJJ - pCfl . The figures are graphs of cross section 
plotted against transverse momentum for various values of COM energy (E cm ) and 
COM scattering angle {0*). Figure |1] shows that the parametrization of Carey et al. is not 
an adequate representation of the data. Figures |] and |3| show that the parametrization of 
Stephens et al. fails for high transverse momentum by severely underpredicting the cross 
section. 

Figure |] shows numerically integrated LIDCS parametrizations of Stephens et al. 
(equation |6|) |pCfl , of Carey et al. (equation U) fll5| , and of equation (|7|) (referred to as 
Kruger) for 7r° production plotted with a parametrization of the integrated formulae of 
Stephens et al. referred to as Stephens-total-param (equation |29[). Three data points 
from Whitmore j24] show that Carey's parametrization does not integrate to the correct 
values and that the rest are quite accurate (see for more detail). 

Figure |^ shows ir° spectral distribution parametrizations given by equations (^) and 
(^|) plotted with LIDCS parametrization of Stephens (equation numerically integrated 
at several lab kinetic energies. Figure |6] shows 7r° spectral distribution parametrizations 
given by equation (|32|) plotted with the numerical integration of equation (0). The shapes 
of the two spectral distributions look quite different even though both original LIDCS 
formula have a similar fit to the data at low p± where the cross section is the greatest, 
and both integrate to the same total cross section. This implies that the available data 
is not sufficient to tightly constrain the shape of the spectral distribution. 

As discussed previously, Figures - |T2] show 7r + and 7r~ LIDCS parametrizations of 
Alper et al. (equation |l"3"D p5[ , Badhwar et al. (equation |16|) ||31|| , Ellis et al. (equation 
|TJ) Igg, Carey et al. (equation |T§ 0, and Mokhov et al. (equation |H]) and LIDCS 
data from [[U], ^| plotted against transverse momentum (p t = p±) for different values of 
COM energy (E cm ), but all at 6* = 90°. These graphs show that the parametrizations 
of Badhwar best fit the data, but underpredict the cross section for large transverse 
momentum. 

Figure [T3| shows the numerically integrated LIDCS parametrizations of Badhwar et 
al. (equation [16]) and of Carey et al. (equation [L4]) |3(J for 7r + and 7r~ plotted 
with parametrizations of the integrated formulae of Badhwar referred to as present work 
(equations |3^ and [31]). Three data points from Whitmore et al. [24) show that Carey's 
parametrization does not integrate to the correct values and that Badhwar's formula is 
accurate. The figures also show that the parametrization fits the numerically integrated 
formulae very well. 

Figures [14] and [15| show 7r~ and rr + spectral distribution parametrizations plotted with 
LIDCS parametrization of Badhwar et al. (equation [16]) |3T| numerically integrated. The 
plot is of cross section (^) plotted against the kinetic energy of the produced pion T n 
at several values for the lab kinetic energies of the colliding proton. The graphs clearly 
show that the spectral distribution parametrizations have excellent fits to the integrated 
LIDCS parametrizations. 
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4. Summary and Conclusions 



This paper presents parametrizations of cross sections for inclusive pion production in 
proton-proton collisions. The cross sections of interest are Lorentz Invariant Differential 
Cross Sections (LIDCS), lab frame spectral distributions, and total cross sections. For 
neutral pions the parametrization of Stephens et. al. [ED] (equation ||) fit the data well 
for low values of p±, but overpredicted the cross section by many orders of magnitude at 
high p± values. Because of this inaccuracy, equation (^) was developed. The final form of 
our resultant parametrization for the neutral pion invariant cross section in proton-proton 
collisions is equation (fjj with D(p±, y/s,9*) given in equation ; F(p±, y/s) given in 



equation (j^), and G(q,p±) given in equation (j^). This formula is as accurate as that 
of Stephens et. al. pDJ at low p± values, but is much more accurate at high p± values. 



For charged pions the formulae of Badhwar et al. (equation [16]) |3T| were found to best 
represent the data except at high p± values. These formulae were used in the development 
of spectral distributions and total cross sections because they are the most accurate at 
low p± where the cross section is the greatest. 

The data for lab frame spectral distributions and total cross sections is scarce, so 
parametrizations for these quantities were developed using the above LIDCS formulae. 
These formulae were numerically integrated, resulting in discrete numerical "data" points 
for these other cross sections. The accuracy of the representations of lab frame spectral 
distributions and total cross sections is, therefore, limited to the accuracy of the original 
LIDCS. The numerical "data" was then parametrized so that closed form expressions 
(equations P3"| - |3"!?D could be obtained. As a check on the accuracy, the total cross section 
numerical "data" was compared to experimental data. They were found to agree quite 
well, but when the numerical "data" for the spectral distributions for the formulae for tt° 
production (equations [23l - p4| and |32]) are compared (ie. compare Figure ^| to Figure |6]), 



they are found to disagree. Since both original LIDCS formulae fit the data well at low 
p± where the cross section is greatest, and both formulae integrate to the correct total 
cross section, the available data must not be sufficient to uniquely determine the global 
behavior of the LIDCS. The data for charged pion production is much more limited than 
the data for neutral pion production, so the same problem exists for charged pions. 

To more accurately determine the cross sections for space radiation applications, mea- 
surements of the spectral distribution at lower energies (for example, proton lab kinetic 
energies 3 and 8 GeV, and pion lab kinetic energies of 0.01 GeV, 0.1 GeV, and 1 GeV) 
would need to be taken. These measurements would put a much tighter constraint on the 
global properties of the LIDCS, and the spectral distribution parametrizations could also 
be made more accurate. 
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Appendix A. Synopsis of Data Transformations 



The data that was used in the comparison of different parametrizations was given in 
terms of several different kinematic variables. Some of the LIDCS data was transformed 
so that all data would be expressed in terms of the same variables: p±, 9*, and E cm = y/s. 
The following is a synopsis of the transformations that were performed for the data plotted 
in the figures. 

The data from Carey et al. 



15] was listed for different values of P p , p±, and 9. \/s, 



p±, and 9 were used by Eggert et al. Stephens et al. [£0| used photon production data 
from Fidecaro et al. [14| to derive pion production cross sections. The variables T} a b, 9, 



and p were used by Stephens, \/s, p±, and the longitudinal rapidity y were used by Alper 
et al. ||25||, but only data with y = was used in the figures. When y = then 9* = 90°. 



The necessary transformations are as follows. 

p± = p sin 9 



(33) 



Ti a b, P p , and 9 can be transformed into y/s and 9* by using the following Lorentz trans- 
formations to change to the COM frame. First express T[ ab and P p as total lab energy 
E. 



E = Ti ab + m p = 
Now perform the following Lorentz transformations. 



PI + ml 



E r 



9* 



jvp cos 9 + jE 
psm 9 



tan 1 



K 7P cos 9 — 'jvE ' 



where 



7 



T tab + 2m r 



7" 



(34) 

(35) 
(36) 

(37) 
(38) 
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P, (GeV) 



Figure 1: 7r° production parametrizations of Carey et al. (equation 5) [|Tf| , of Stephens et 
al. (equation 6) |2(J, and of equation (7) plotted with LIDCS data from |§ gOj. LIDCS 
is plotted against transverse momentum for COM energy E cm = 7 GeV. The pion COM 
scattering angle is 12.2 < 6* < 12.4° and 6* = 89° for the data, and the parametrizations 
are plotted at 6* = 12° and 6* = 89°. 
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Figure 2: 7r° production parametrizations of Carey et al. (equation 5) of Stephens et 
al. (equation 6) p0| , and of equation (7) plotted with LIDCS data from |L3|, 0, [17|| . The 
data is for COM energy E cm = 62.4 GeV and 62.9 GeV at COM scattering angle 9* = 90° 
and the parametrizations are plotted at E cm = 62.6 GeV at 9* = 90°. The second set of 
data is at COM energy E cm = 23.3 GeV and the pion COM scattering angle 9* = 15° 
and 17.5° for the data, and the parametrizations are plotted at 9* = 16°. 
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Figure 3: n° production parametrizations of Carey et al. (equation 5) jl5| , of Stephens et 
al. (equation 6) p0| , and of equation (7) plotted with LIDCS data from ([HJ. LIDCS is 
plotted against transverse momentum for COM energy E cm = 18.2 GeV. The pion COM 
scattering angle 32.3° < 9* < 32.5° for the data, and the parametrizations are plotted at 
6* = 32°. The second set of data and parametrizations are at E cm = 53 GeV, 6* = 90°. 
The data is from 113, |16|, 0, ITBl. 
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Figure 4: Parametrization of total n° production cross section plotted with numerically 
integrated LIDCS parametrizations of Stephens et al. (equation 6) [p0[] , of Carey et al. 
(equation 5) fll5l , and of equation(7) referred to as Kruger. The curve labelled 'Stephens- 
total-param' is the parametrization given in equation (29). Three data points from Whit- 
more 



24] are included for comparison. 



50 




T„ (GeV) 

Figure 5: 7r° spectral distribution parametrizations of equations (23) and (24) (solid lines) 
plotted with LIDCS parametrization of Stephens (equation 6) pO) numerically integrated 
at lab kinetic energies of 0.5 GeV, 1.0 GeV, 1.9 GeV, 5.0 GeV, 9.5 GeV, 20 GeV, and 50 
GeV, listed in order of increasing cross section (symbols). 



26 



60 




T„ (GeV) 

Figure 6: n° spectral distribution parametrization of equation (32) (solid lines) plotted 
with equation(7) numerically integrated at lab kinetic energies of 0.5 GeV, 1.0 GeV, 1.9 
GeV, 5.0 GeV, 9.5 GeV, 20 GeV, and 50 GeV, listed in order of increasing cross section 
(symbols). 
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Figure 7: n~ production parametrizations of Ellis et al. (equation 15) p2[ , Badhwar et 
al. (equation 16) Carey et al. (equation 14) |HJ, and of Mokhov et al. (equation 17) 
|32|1 plotted with LIDCS data from LIDCS is plotted against transverse momentum 
for COM energy E cm = 23 GeV and pion COM scattering angle 9* = 90° 
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Figure 8: n production. Same as Figure 7 except E cm = 31 GeV, and the parametrization 
of Alper et al. (equation 13) is included. 
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Figure 9: 7r + production parametrizations of Ellis et al. (equation 15) p2] , Badhwar et 
al. (equation 16) |^TJ and of Mokhov et al. (equation 17) |32| plotted with LIDCS data 
from p5| . LIDCS is plotted against transverse momentum for COM energy E cm = 23 
GeV and pion COM scattering angle 9* = 90° 
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Figure 10: ix + production. Same as Figure 9 except E cm = 31 GeV, and the parametriza- 
tion of Alper et al. (equation 13) p5f is included. 



10" 



CD 

CD 



1 O 



1 O 



1 O 



1 O 



* E cm =45 GeV 0=9O° 

Badhwar 

Mokhov 



10° 

P t (GeV) 



\ 

\ 



10' 



Figure 11: n + production. Same as Figure 9 except the data is from (T^, ^5) for E cm = 
45.0 GeV and 44.8 GeV, and the parametrization of Ellis et al. (equation 15) [E| is 
excluded. Parametrizations are plotted at E cm = 45.0 GeV. 
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Figure 12: 7r~ production. Same as Figure 7 except the data is from |Tc], go) for E cm = 
45.0 GeV and 44.8 GeV, and some of the parametrizations are excluded. Parametrizations 
are plotted at E cm = 45.0 GeV. 
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Figure 13: Parametrizations of total 7r =l= production cross section (present work) (equations 
30 and 31) plotted with numerically integrated LIDCS parametrizations of Badhwar et 
al. (equation 16) |3l[] (circles and squares) and Carey et al. (equation 14) |30| (triangles). 
Six data points are included for comparison (data is from Whitmore WW ). 
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Figure 14: ir~ spectral distribution parametrizations (equations 27 and 28) (solid lines) 
plotted with LIDCS parametrization of Badhwar et al. (equation 16) |3l| (symbols) 
numerically integrated at lab kinetic energies of 0.5 GeV, 1.9 GeV, 5.0 GeV, 9.5 GeV, 20 
GeV, and 50 GeV, listed in order of increasing cross section. 
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Figure 15: tt + spectral distribution parametrizations (equations 25 and 26) (solid lines) 
plotted with LIDCS parametrization of Badhwar et al. (equation 16) |JT| (symbols) 
numerically integrated at lab kinetic energies of 0.5 GeV, 1.1 GeV, 5.0 GeV, and 50 GeV, 
listed in order of increasing cross section. 
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